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Abstract 

Most analyses of dark matter within supersymmetry assume the entire cold dark matter 
arising only from weakly interacting neutralinos. We study a new class of models consisting of 
?7(1)" hidden sector extensions of the minimal supersymmetric standard model that includes 
several stable particles, both fermionic and bosonic, which can be interpreted as constituents 
of dark matter. In one such class of models, dark matter is made up of both a Majorana dark 
matter particle, i.e., a neutralino, and a Dirac fermion with the current relic density of dark 
matter as given by WMAP being composed of the relic density of the two species. These 
models can explain the PAMELA positron data and are consistent with the antiproton flux 
data, as well as the photon data from FERMI-LAT. Further, it is shown that such models 
can also simultaneously produce spin-independent cross sections which can be probed in 
CDMS-II, XENON-100 and other ongoing dark matter experiments. The implications of the 
models at the LHC and at the next linear collider (NLC) are also briefly discussed. 
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1 Introduction 



Recently several particle physics models have been constructed that connect the standard 
model (SM) to hidden sectors and lead to massive narrow vector boson resonances as well 
as other signatures which can be detected at colliders [H 121 [3]. The connection to the 
hidden sector arises via mass mixings and kinetic mixings [H [21 El IH El E] and via higher 
dimensional operators. Models with the above forms of communication between the sectors 
also have important implications for dark matter [71 [3l [6] (for a review see [El |9]). In 
this work we show that multicomponent dark matter can arise from f/(l)" extensions of 
the minimal supersymmetric standard model (MSSM) with Abelian hidden sectors which 
include hidden sector matter. Our motivation stems in part from the results of several dark 
matter experiments that have recently appeared. Thus the PAMELA Collaboration [10] has 
observed a positron excess improving previous results from HEAT and AMS experiments 
[TT] . One possible explanation of such an excess is via the annihilation of dark matter in the 
galaxy [T2]. Additionally, recent data from CDMS-II hints at the possibility of dark matter 
events above the background, and this will be explored further by the upgraded XENON 
experiment [131 E] • 

For a thermal relic, the PAMELA data and CDMS-II data taken together at face value 
do raise a theoretical puzzle if indeed both signals arise from the annihilation of cold dark 
matter. Thus most models which aim to explain the PAMELA positron excess do not give 
a significant number of dark matter events in the direct detection experiments currently 
operating. Conversely, models which can give a detectable signal in direct detection exper- 
iments typically do not explain the PAMELA data without the use of enormous so-called 
boost factors. As we will show here, this can be circumvented in models where the dark 
matter has several components. Thus, motivated in part by the recent cosmic anomalies we 
develop supersymmetric models which contain minimally a hidden Abelian sector broken at 
the sub-TeV scale where the mass generation of the hidden states involves nontrivial mixings 
with the field content of the electroweak sector of the minimal supersymmetric extension of 
the standard model leading to dark matter which can have several components which can 
be both bosonic and fermionic. 

More specifically, in this work we go beyond the simple theoretical construction that ther- 
mal dark matter compatible with WMAP observations is composed of a single fundamental 
particle. There is no overriding principle that requires such a restriction, and nonbaryonic 
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dark matter (DM) may indeed be constituted of several components, so in general one has 
{yth?')£,M = '^i{,^h'^)DMi, where i refers to the various species of dark particles that can 
contribute to the total nonbaryonic {Qh'^)DM- In fact we already know that neutrinos do 
contribute to dark matter although their contribution is relatively small. Thus we propose 
here a new class of multicomponent cold dark matter models in Abelian f/(l) extensions of 
MSSM which can simultaneously provide an explanation of the PAMELA and WMAP data 
through a Breit-Wigner enhancement [12], while producing detectable signals for the direct 
searches for dark matter with CDMS/XENON and other dark matter experiments. 

A simultaneous satisfaction of the PAMELA positron excess and the satisfaction of 
WMAP relic density constraints can also occur if there is a nonthermal mechanism for 
the annihilation of dark matter with a wino lightest (R parity odd) supersymmetric particle 
(LSP) [ISl [m [HI HSl El Ej • However, a detectable spin-independent cross section in such a 
nonthermal framework does require that a pure wino is supplemented by a suitable admix- 
ture of Higgsino content as in the analysis of [19] and in [20], the later for a thermal relic. 
We remark that multiple U{1) factors and its influence on dark matter have very recently 
been studied [20l [2T] . We also remark, some other works have recently looked at dark matter 
with more than 1 component [22]. The models proposed and analyzed here are very different 
from these. 

The outline of the rest of the paper is as follows: In Sec. ([2]) we give a detailed description 
of the two models one of which is based on a U{l)x extension of the MSSM where U{l)x 
is a hidden sector gauge group with Dirac fermions in the hidden sector. This model allows 
for dark matter consisting of Dirac, Majorana, and spin zero particles. The second model is 
based on a f/(l)x x U (l)c extension of MSSM, where U (l)c is a gauged leptophilic symmetry 
and U{l)x, as before, is the hidden sector gauge group which also contains Dirac particles in 
the hidden sector. This model too has Dirac, Majorana, and spin zero particles as possible 
dark matter. In both cases we will primarily focus on the possibility that dark matter 
consists of Dirac and Majorana particles, and we will not discuss in detail the possibility 
of dark matter with bosonic degrees of freedom. In Sec.(|3]) we discuss the relic densities in 
the two component models. In Sec. (HI) we give an analysis of the positron, antiproton, and 
photon fluxes in the two models. In Sec. ([5]) we give an analysis of event rates for the proposed 
models for CDMS-II and for XENON-100. We give the analysis within the framework of 
supergravity grand unified models [231 |2l] defined by the parameters mo, mi/2, Aq, tan/3, and 
sign(yu) with nonuniversalities (NUSUGRA) defined by 61^2,3 in the gaugino sector so that 
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U{1)y X SU{2)l X SU{?>)c gaugino masses at the grand unified theory (GUT) scale are given 
by fhi = 777,1/2(1 + Si) {i = 1, 2, 3) (see, e.g., [25] and references therein). We also discuss the 
possible new physics one might observe at the LHC (for a recent review see also [9]) and 
elsewhere for these models. Conclusions are given in Seed?]). 

2 Multicomponent Hidden Sector Models 
2.1 Multicomponent U{l)x model 

A U{l)x extension of the minimal supersymmetric standard model involves the coupling of 
a Stueckelberg chiral multiplet S = {p + ia, xs, -^s) to vector supermultiplets X, B, where 
p is a real scalar and a is an axionic pseudoscalar. Here X is the U{l)x vector multiplet 
which is neutral with respect to the SM gauge group with components X = {X^, Xx, Dx), 
and B is the f/(l)y vector multiplet with components {B^, \b, Db), where the components 
are written in the Wess-Zumino gauge. The chiral multiplet S transforms under both U{l)x 
and U{1)y and acts as the connector sector between the visible and the hidden sectors. The 
total Lagrangian of the system is given by 

^ = -^^MSSM + ^Uil)x + -^st (1) 

where is the kinetic energy piece for the X vector multiplet and Cst is the supersym- 

metric Stueckelberg mixing between the X and the B vector multiplets so that [H [7] (see 
also [261 EH ED]) 

Cst = j d^0£9 (MiX + M2B + S + SY , (2) 

where Mi and M2 are mass parameters. The Lagrangian of Eq.([T]) is invariant under the 
U{1)y and U{l)x gauge transformations, i.e., under 

SxX = Cx + Cx, SxS=-MiCx, 6yB = Cy + (y, SyS = -M^Cy, (3) 

where ( is an infinitesimal transformation chiral superfield. In component form we have for 
the Stueckelberg sector with U{l)x x U{1)y 

Cst = -^{MiX, + M2B, + d,af-^{d^pf-ixscr^d,xs + 2\Fs\' 

+p{MiDx + M2DB) + Xs{Mi\x + M2\b) + Xs{MiXx + MsAs) . (4) 
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In addition, one may include a supersymmetric kinetic mixing term between the U{l)x and 
f/(l)y gauge fields [7] leading to £ = £mssm + £(7(i)x + -^km + £st, where 

-^X'^'^B^, - tdiXxa^'d^XB + Xscr'^d^Xx) + 6DbDx • (5) 

One can also add additional D terms as in [7j. Both Stueckelberg and kinetic mixings of the 
gauge fields U{l)x and f/(l)y are constrained by the electroweak data[2]. As a consequence 
of the mixings, the extra gauge boson of the hidden sector couples with the standard model 
fermions and can become visible at colliders. The Lagrangian for matter interacting with 
the U{1) gauge fields is given by 

£matt = / dW J2 [<|.,e2^-«-^+2'^-«-^$, + •Ihid.e^'^-Q-^+^^-Q-^Khid.] . (6) 

where the visible sector chiral superfields are denoted by $i (quarks, squarks, leptons, slep- 
tons, Higgs, and Higgsinos of the MSSM) and the hidden sector chiral superfields are denoted 
by $hid,j- In the above, Qy is the hypercharge normalized so that Q = T^+Qy- As mentioned 
already, the SM matter fields do not carry any charge under the hidden gauge group and vice 
versa, i.e. Qx^i = and QsM^wd = 0. The minimal matter content of the hidden sector 
consists of a left chiral multiplet $hid = (0)/)-^) and a charge conjugate $1;;^ = 
so that $hid and ^^^^ have opposite U{l)x charges and form an anomaly-free combina- 
tion. A mass for the Dirac field tp arises from an additional term in the superpotential 
= M^^^'^, where ip is composed of / and /'. The scalar fields acquire soft masses of 
size rriQ from spontaneous breaking of supersymmetry by gravity mediation, and in addition 
acquire a mass from the term in the superpotential so that 

ml = ml + Ml = ml,. (7) 

After spontaneous breaking of the electroweak symmetry there would be mixing between 
the vector fields X^, B^, A^fj^, where A^^ is the third component of the SU{2)l field Aa^, 
(a = 1,2,3). After diagonalization V'^ = {X,B,A3) can be expressed in the terms of the 
mass eigenstates E'^ = {Z', Z, 7) as follows: 

V^ = 0,jE,, i,j = 1-3, E={Z',Z,^). (8) 
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The neutral vector mass squared matrix is of the form given in Ref. [1] of [B]. Further, the 
chiral fermions in the S + S multiplet together with the MSSM gauginos and Higgsinos will 
form a 6 X 6 neutralino mass matrix whose eigenstates are six neutralino states Xa, cl = 1 — 6, 
where we assume that the set Xi ■ ■ - xt is the regular set of neutralinos and x% Xe ^^^^ the 
two additional neutralinos that arise in the U{l)x extension. From the components Ax, Ax 
and xs, Xs that appear in Eq.(j4]), we can form two Majorana fields Ax and ips as follows: 



These components combine with the MSSM gauginos and Higgsinos to form a 6 x 6 neutralino 
mass matrix whose eigenstates are the six neutralinos Xa, (a = 1 — 6). Thus Ax and ips can 
be expanded as linear combination of Xa, i-e.. 



where R is the unitary matrix which diagonalizes the 6x6 neutralino mass matrix. Further 
the CP even Higgs sector is extended by the additional state p [T] . The results outlined here 
give the following types of interactions: 

1. There are interactions of the Dirac fermion in the hidden sector with the standard 
model particles via Z,Z', 7 interactions. Thus, the Dirac dark matter can annihilate 
into standard model particles via exchange of Z,Z', 7 in the early universe and in the 
galaxy. Depending on which of the two, Dirac or Majorana, is the heavier one may 
have Dirac particles annihilating into Majoranas or the Majorana particles annihilating 
into Dirac fermions in the galaxy: 



2. In addition to the above we have fermion-neutralino-sfermion couplings in the hidden 
sector as given by Eq.(l6]). Thus interactions of the type ipXa<P + h.c, etc. can produce 
decays such as — )■ + Xa if they are kinematically allowed. 

3. The scalar field p is CP even and mixes with the MSSM Higgs fields. Through these 
mixings p has couplings to the SM fermions and through these couplings it can decay 
into the SM fermions. 




(9) 



Ax = RlaXa, a = 1 - 6, Ips = R2aXa, = 1-6 



(10) 



i^iJ XX or XX • 



(11) 
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It is instructive to list all the new particles in this U{l)x model as summarized below: 

New particles of the U{l)x model 
spin : p, 0, 0', 

spin ^ : i^,xl,xl, 

spin 1 : Z'. (12) 

We assume that the lightest i?-parity odd particle (LSP) is the least massive neutralino 
(x° = Xi = x) cind resides in the visible sector and thus the masses of xt^ Xe ^ire larger than 
the LSP x° mass, and consequently Xs, Xe are unstable and decay into SM particles and 
The bosons Z' and p are unstable and decay into SM fermion pairs // with the decay of the 
p going dominantly through the process p bb or p ^ tt ii nip > 2mt. The remaining three 
particles 0, 0' are all milli charged and, consequently, at least one of them is stable. If 
we assume m<^, m^i > M^, at least ip is always stable and the other two may or may not be 
stable. These along with the LSP give rise to various possible candidates for dark matter. 
Thus, depending on the relative masses of the Majorana, Dirac, and spin particles there 
are three possibilities for the constituents of dark matter as outlined below. 

Two component dark matter: Majorana + Dirac. This model arises as follows: 
consider the case where > + M^. In this case the decays 0,0' — )■ + x"? will 
occur and 0, 0' will be unstable. Thus ip is stable and so is x under the assumption of R 
parity conservation. Consequently, we will have two dark matter particles; namely, one a 
Majorana which is the LSP in the visible sector and the other a Dirac in the hidden sector. 
The Majorana and Dirac particles once created will annihilate as follows: 

ij + ^jj ^ Z,Z',-f ^SM + SW, (13) 

X + X^{s: Z', Z, h, H, A, p), {t/u : xu Xk) ^ SM + SM'. (14) 

where s : and t/u : refer to s and t or u channel exchanges. In addition to Eq.( |T4l) there 
are coannihilation processes which contribute to the relic density. Since both ip and x ^i-re 
stable, the total relic density of dark matter will be the sum of the relic densities for the two, 
the sum being constrained by the WMAP data. These constraints are discussed further in 
SecilaD. 



7 



Three component dark matter: Dirac and two spin particles. Suppose the mass 
of X is larger than the sum of the masses of the Dirac plus the scalar 0, i.e., M-^ > M^ + m^. 
In this case the decay x — ^ + V^, 0' + will occur and, consequently, x is unstable. On the 
other hand, 0, 0' and ip are stable since they cannot decay into anything else. Thus, here we 
have three dark matter particles: one Dirac, and the other two spin 0. Processes that lead 
to the annihilations of these particles are those in Eq. (fT3|) for ip, and also for and 0', they 
are similar to those in Eq.([T3D, i.e., + 0*, 0' + 0'* 7, Z, Z' ^ SM + SM'. In this three 
component dark matter model all the components reside in the hidden sector and thus their 
couplings to the standard model particles are extra weak. Consequently, they will have very 
small spin-independent cross sections in direct detection experiments. For this reason, this 
class of models is less preferred compared to the two component model. 

Four component dark matter: Majorana, Dirac, and two spin particles. Finally, 
we consider the case when either of the following two situations occur: (i) > M^jfri^ < 
< + m^, (ii) < < + M^. In these cases all four particles, one Majorana, 
one Dirac, and two spin particles, are stable and thus are possible dark matter candidates. 
These particles will annihilate to the SM particles as in Eq. (fT3|) . Eq.f ll4p and for and 0' 
via processes in the three component dark matter model as described above. This model is 
in many ways similar to the two component model and like the two component model this 
model too should lead to detectable signals in experiments for the direct detection of dark 
matter. 

2.2 Multicomponent Leptophilic U{l)x x U{l)c model 

We discuss now another model which contains two additional Abelian vector bosons where 
one of the extra bosons is leptophilic. Leptophilic Z's have a long history [28] and have been 
revisited [22] over the recent past in the context of dark matter. Here we will consider a 
U{l)x X U{l)c model where the U{l)x as before is in the hidden sector, and U{l)c is a 
leptophilic symmetry. As in the U{l)x model, we also assume that the hidden sector has 
a pair of Dirac fermions ip and -0 which are charged under U{l)x but are neutral under 
the standard model gauge group and under U{l)c- Regarding f/(l)c "we assume it to be 
Lg — L^, i.e., a difference of family-lepton numbers, which is anomaly free, and can be gauged. 
The corresponding gauge field couples only to e, fj, families and nothing else. The total 
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Lagrangian in this case is 

1^ = -^^MSSM + ^U{1)'^ + -^St, (15) 

where Cu{i)'^ is the kinetic energy for the X and C multiplets and for £st we assume the 
following form: 



d^ed^e {MiC + M'^X + M'^B + s + sf 

+ j d^ed^e iM[c + M2X + m;;b + s' + s'f 



(16) 



where C is the U{1)l^-l^ vector multiplet with components {C^, \c, Dc) and X and B 
are the U{l)x and U{1)y multiplets as discussed before. The gauge transformations under 
U{l)c, f/(l)x, and U{1)y are 



ScC = Cc + Cc, ScS=-M,Cc, ScS' 
SxX = (x + Cx , 5xS = -M'^Cx , 5xS' = 
6yB = Cy + Cy , SyS = -M'^Cy , ^yS' = 



-M[Cc 
-M2CX , 



(17) 



where (c, Cx, Cy, etc. are the infinitesimal transformation chiral superfields. The quantities 
Ml, M2, M(, MgjMg, and M3 are the mass parameters. In the vector boson sector £st 
assumes the form 



Cst = --{M,C, + M^X^ 



The mass^ matrix in the vector boson sector in the basis (C^, X^, 5^, A'^^) is given by 



/ M^ + Mf MiM!^ + M[M2 MiM'^ + M[M'^ 

MiM!^ + M[M2 M| + M!^ M^M^ + MaM^' 

MiM;^ + M(M^' M!^M'^ + M2M'^ M^^ + + 

\ -MyMw 






-MyMw 



(19) 



where My[r=g2-v/2 is the W boson mass and My = Myy ta.n 6y[r= gy ■ v/2, and where 
6w is the weak angle. The dynamics of the model of Eq. f|T9|) is rather involved. We will 
focus, therefore, on a simpler version of this more general case where we neglect the mixings 
with B^, i.e., we set M3 = M3 = 0. Inclusion of these coupling in the analysis would not 
drastically change the analysis or the conclusions of this work as long as we keep the mixing 
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parameters M3/M1 2, Mg /Mi 2 very small. After neglecting the mixings with 5^, the mass^ 
matrix is block diagonal and so we can diagonalize the top left hand corner 2x2 mass matrix 
independent of the standard model sector. We are interested in the limit of small mixing 
between U{l)x and U{l)c and thus consider^ 



M[,M^ < Mi,M2. 
In the above approximation the eigenvalues of this mass matrix are 

Am2 ~ 



(20) 



(MiM^ + M[M2f 



[Ml + - Ml - M'iY 
The corresponding mass eigenstates are Z' and Z" , where 



(21) 



= cosOxZ'' — sm6xZ' = sinOxZ'' + cos OxZ' 



tan 9 



X 



MiM^ + M[M2 



(22) 



Because of Eq. fl20p tan^x ^ 1- In the above, the Dirac fermions in the hidden sector have 
no couplings with the photon and are electrically neutral. However, by a small mixing of 

with in Eq. (fT8|) . we can generate a milli charge for the Dirac particles in the hidden 
sector consistent with all electroweak data. 

We discuss now the gaugino/chiral fermions in the extra f/(l) sectors which arise from 
the superfields C,X, S + S, S' + 5". From the gaugino components Ac, Ac, Ax, Ax, and from 
the chiral fermion components in the extra U{1) sectors Xs^Xs^Xs'^Xs'^ one can construct 
four component Majorana spinors two of which are exhibited in Eq.([9]) and the remaining 
two are given by 



A, 



c 



Aca 



Xa,S' 



(23) 



The neutralino mass matrix in the [U{l)x x f/(l)c] x [SU{?>)c x SU{2)l x U{X)y] model 
takes a block diagonal form 



" f/(l)x X f/(l)c 


04x4 


sector 




04x4 


MSSM 




sector 



(24) 



^x8 



^Note these mass terms il/i,M2 are different than those considered in Sec. 12.1 
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Thus, the Stueckelberg mass generation produces a mass matrix in the hidden gaugino/chiral 
fermion sector which is decoupled from the neutralino mass matrix in the visible sector. 
Specifically in the 4 component notation the gaugino/chiral fermion mass matrix in the 
U{l)x X U{l)c sector is given by 



T mass 

^U{l)xxU{l)c 







T 


/ 








Ml 


K \ 






\ 
















M[ 










Ac 








Ml 


M[ 










Ac 




\ Ax 


J 




v 




M2 





J 




V Ax 


/ 



(25) 



In the diagonalized basis we can label the extra neutralinos by Xg, Xg, x?, Xs- Since the 
hidden sector and the neutralinos of the visible sector are decoupled, the diagonalization of 
the neutralinos in the visible sector, i.e., of x^, = 1 — 4) is not affected. Further, as for the 
case of the U{l)x model, it is instructive to list all the new particles in this U{l)x x U{l)c 
model as summarized below: 

New particles of U{l)c x f/(l)xniodel 
spin : p^p^cp^cj)', 

spin - : ^A, Xs' Xe, X?, xl. 

spin 1 : Z\ Z" . (26) 

We discuss now the stability of the new particles in this model. As before we assume that 
the mass of (and of 0') is larger than the mass of ip. Thus ip will be stable since it cannot 
decay into anything. If kinematically allowed the fields (p and (p' can decay only via the 
process (p,(p' ^ ip + ^is in the U{l)x model. Of the remaining fields obviously Z' and Z" 
are unstable as they decay into ee,pfi, u^Ue, u^u^ as well as into -^-^ depending on the mass 
oi ip. As already noted, a small milli charge can develop for the hidden sector matter via 
small couplings of the and fields. The phenomenology of such models will be very 
similar to the one we are discussing here. 

The extra neutralinos of Eq. fl26|) can also be all unstable. Thus h^c couples with leptons- 
sleptons (e, e etc.) via coupling of the type AcCie^ , etc. and after diagonalization of the 
gaugino/chiral fermion mass matrix all the x^? (/c = 5 — 8) will have coupling with leptons- 
sleptons of the type indicated. Further, two of the Xfc have roughly a mass of size Mi while 
the remaining two have roughly a mass of size M2. Thus, if Mi,M2 > m-yA), which is what 
is assumed in this work, all the neutralinos of the hidden sector will be unstable and decay 
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into final states of the type eex°,/i/ux°, etc. Regarding the field p, there is an interaction of 
type 



With this interaction p will decay as follows: p — )■ /*/ — )■ ffx^X^if = provided this 
process is kinematically allowed which we assume is the case. A similar situation occurs for 
the case of p'. Additionally, if there is a mixing with in the Stueckelberg sector then, as 
in the analysis of the U{l)x model, the fields p and p' will mix with the Higgs sector and can 
have decays of the type p bb, p' — t- bb, etc . Thus, in the end we are left with a similar 
set of possibilities for dark matter as in the U{l)x model, i.e., (i) a two component model 
with ip and x^j (h) a three component model with ip, 0, 0', and (iii) a four component model 
with 0, 0', and However, as in the U{l)x case we will focus on the two component 
model consisting of Dirac and Majorana dark particles. 

We assume M|„3>M|, and that the annihilation of dark matter occurs close to the Z' 
pole for reasons that will become apparent shortly. As a consequence, the annihilation of 
dark matter in the early universe and in the galaxy is controlled by the Z' pole and the effect 
of the Z" pole on the analysis is essentially negligible. The basic interaction of and of 
with matter is given by 



where / runs over e and p families and where Qq = —Qq- the mass diagonal basis the 
interaction of Eq. fl28|) assumes the form 



The interaction of Eq.f l29|) leads to the annihilation of ipil) into e^e~ and p~^fi~ via the Z', Z" 
poles for which we assume a Breit-Wigner form. Thus, the ■0-0 — ?■ // annihilation cross 
section takes the form 



M,gcpirQy~), f = e,p. 



(27) 



(28) 



Ant = (gxQxi^l^ip COS 9 X - gcQcfl^f sm9x)Z'^ 

+{gxQxhy sin Ox + gcQUrf cos ex)z';^. 



(29) 



a^\{s-Ml, + iTz'Mz')- 
(3f{gxgcQxQ{,sm{29x)) 

64:713(3^ 



-1 
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where = (1 — Amj ^/ sY^'^. The relevant partial Z' decay widths are given by 

nZ'^ff) = {gcQf^smex)'^, / = e,/x, (32) 

i/vr 



Y{Z' ^ = [g^Qx^o^Q^f-^ [l + [\ - e(Mz. - 2M^), (33) 

and similarly for the partial decay widths of the Z" with M^/ — )■ Mz» and — sin 6*^ — )■ cos 6*^ 
in Eq.([32D and cos^x ^ sin^x in Eq.dSS]). 

A constraint on gc comes from the contribution of the Z' and Z" to g^ — 2 [5U]. Their 
exchange gives 



247r2 



sin^ 6x cos^ 6'x 



(34) 



Using the current error [30] of A((7^ — 2) = 1.2 x 10 ^ in the determination of g^ — 2 and 
assuming 6'x is small, one finds the following constraint on ac'- 

ac < 0.001 J — , 35 

where ac = g'c/^'^- We note that if the mixing angle dx is small, the decay width of Z' — )■ // 
(/ = e,/i) and of Z" — ipip will be narrow while the decay width of Z" — t- // (/ = e,/i) 
and of Z' — )■ ipil} will be of normal size. However, when ~ M^//2 the Z' decay width 
into '?/^'?/^ will also be small due to the kinematic suppression factor (1- [AMU Ml]) ' . 
this case we will have the total width of the Z' to be rather narrow. Thus for annihilation 
near the Breit-Wigner pole we will have a large enhancement of {av) due to the narrowness 
of the Z' [12]. It was shown in the analysis of Feldman-Liu-Nath in [12] that near the 
Breit-Wigner pole such annihilations allow one to fit the relic density as well as allow an 
enhancement of {av) in the galaxy. We note that while Z' decay width is very small this is not 
necessarily the case for Z" which can decay into ee, t'e'^e? v^J.'^^^l with normal strength. Thus 
neglecting the contribution of Z" — ipip which is small due to the sin^ 6x ~ suppression, 
one finds the total width of Z" to be T z" — co^'^ 9 xotcMz"- We will see in Sec.(jlj) that 
the ac needed in the analysis of the relic density is relatively small compared to normal 
electroweak coupling and, consequently, the width of Z" though significantly larger than the 
Z' width is still relatively small compared to what one might expect for a Z' in a GUT model 
and certainly much smaller than the width for a Z' arising as a Kaluza-Klein excitation in 
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the compactification of an extra dimension [HH [32] • Finally, the annihilation of the Dirac 
particles in the early universe goes by the processes 

ipip Z', Z" -> e^e~,fi^fi', UeUe, '^ii'^ii, (36) 

which is to be contrasted with the processes Eq. (fT3l) in the U{l)x model. 



3 Relic Density in a Two Component Model 

Here we discuss the relic density in models with two components. A general analysis requires 
solving the Boltzmann equations in a Friedmann- Robert son- Walker universe [331 IM], and 
includes coannihilations [32] and an accurate integration over pole regions. As in the MSSM 
alone, one will generally encounter the Z and Higgs poles [M] and these need to be treated 
with care. The number changing processes include 

ipij O SM SM', -^ij o XX, XX ^ SM SM'. (37) 

Note that the process ipx ^ SM SM' is not allowed since ipx connect only to and 0', 
neither of which can connect to the standard model particles. For the simplest two component 
model with dark matter particles ip, x, with the assumption that > the only relevant 
processes in the annihilation of ipip are ipip — )■ //, XX final states. Since ip is heavier than 
X its freeze-out occurs earlier (at a higher T) than for x- Thus, the Boltzmann equations 
for (which includes fermions and antifermions) and for for the U{l)x and for the 
U{l)x X U{l)c two component models are given by 



-3Hn^- -{av)^^{nl-nl^^^), (38) 



— - -SHn^ - {(Tv)^^{n^ - n^^J + -{crv)^^_^^^{n^ - n^,q). (39) 

Here (av)^^ refers to ipip //, XX, and {av)^,^ stands for {av)^^^^^^ For the spin 
averaged cross section for the Dirac case, the extra factor of 1/2 is to account for the fact 
that we are dealing with a Dirac fermion. The number densities are n^, and n^^eq, ^.eq are 
their values at equilibrium, i.e., n(^^^^)^cq — fl'(i/;,x)(^(i/',x)-^)/2^)^^^^^P(~~~T^)' "W'here g^p = 4: 
and = 2. Since the two dark matter particles are sub-TeV in mass, they will freeze-out at 
temperatures that are not drastically different. One can solve the Boltzmann equation for ip 
with the appropriate boundary conditions to compute the freeze-out temperature and the 
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relic density of ip at the current temperatures. To compute the freeze-out temperature Tj 
for the particles Xi o^is uses solutions for as computed from the Boltzmann equation for 
ip as input in the Boltzmann equation for x keeping in mind that n^,eq in the x Boltzmann 
equation can be neglected since we are below the freeze-out temperature for ip. It is difficult 
to get a closed form solution of Eq. (15^]) for and thus in general the analysis must be done 
numerically for Vt^h?. However, it turns out that for both the U{l)x and the U{l)x x f/(l)c 
models the contribution of the term proportional to in Eq.f l39l) is rather suppressed and 
it is a good approximation to neglect this term for both models. In this case, one has 

(l^/i')wMAP = {n^h\ + {n^h\ ^ % + %, (40) 



where 



1.07 X 10^ GeV-i ^ 1.07 X 10^ GeV-i 
~ ; , Q ~ 2 X , , (41) 



Jo = / (^^').. dx , 4= (at;) , dx , (42) 



Jo 

and where g*{ip, x) denotes the effective degrees of freedom at the freeze-out of tp, x respec- 
tively. The analysis leading to Eqs. ( 138|39|) is easily extended to include coannihilations. The 
analysis can easily be reversed if the Majorana is heavier than the Dirac. Denoting p©,^, 
Pq^-^ as the local density of each dark matter kind in the halo, one can assume 

P0,v/P0,x ~ (l^^/i^)o/(^^x^^)o- (43) 

However, the ratios need not be the same. The local halo densities are also constrained such 
that Pq^^ + Pq^^ = P0,totai — (0.35 — 0.45)GeVcm~^. For the calculation near the Z' pole, we 
use the analysis of of [6] which follows the techniques of [36] . Indeed the analytic techniques 
developed in [361 [6] have been cross-checked with independent codes. For the U{l)x model, 
the decay branching ratios are substantially less hadronic and more leptonic than for the 
annihilations via the Z boson exchange [37]. For the U{l)x x U{l)c model the decays of the 
Z', Z" are purely leptonic. These leptophilic decay patterns for the extra Z's help to explain 
the PAMELA positron excess without recourse to large ad hoc boost factors. 
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4 Positron, Antiproton, and Photon fluxes in the U{l)x 

and U{l)x x U{l)c Models 

An excess of positrons, antiprotons and photons over the cosmic background is a possible 
indicator for annihilating dark matter in the galaxy as was pointed out early on [38] . In our 
multicomponent supersymmetric models contributions from the Majorana component to the 
fluxes are negligible (suppressed by an order of magnitude or more) and essentially the entire 
effect arises from the Dirac component. The positron flux $£+ arising from the annihilation 
of dark matter (DM) particles is [391 HQ] 

Here M^m is the mass of the dark matter particle. In the above rj = (1/2, 1/4) for (Majorana, 
Dirac) cases, respectively, Bg is a boost factor which may arise as a consequence of dark 
matter substructure, or a local clump. Recent N-body simulations show that it is unlikely 
that large dark matter clumps exist within the halo of our galaxy [H] and thus the use of 
large clump factors in flux analyses appear unreasonable [12]. 

The other parameters that enter Eq.( l44l) are as follows: X(^e,e') is the halo function and we 
parametrize it in some of the standard forms adopted in the literature with the appropriate 
diffusion models [39] using the standard profiles [13]. The positron velocity is Vg.+ ~ c, and 
the energy loss function b{E) has the form b{E) = Eo{E / EqY /te-, where Te ~ (1 — 2)10^^[s], 
with E in [GeV] and Eq = IGeV. We use the GALPROP background estimate of [44] fit in 
[l5] ; modifications of the background estimates require either smaller or larger mass splitting 
at the pole which can range from the order a GeV to the order of tens of MeV, depending on 
the assumed astrophysical background and the level of dumpiness of the signal. Further, 
in Eq. p^ ('^^)haio velocity averaged cross section in the halo of the galaxy. We note 

that (o''y)haio significantly different than {<^v),^^ at the epoch of freeze-out. Thus as 

emphasized in [12] a replacement of {o'v)^^^^ by freeze-out' often done, is inaccurate 
and can lead to significant errors in the positron flux computation. This stems from the fact 
that the relic density in previous works is often approximated by pulling out (erf) from the 
integral between Tfreezc-out and the current temperature. In general, full integration must be 
taken into account in the vicinity of a pole for an accurate calculation [36] or when the dark 
matter coannihilates [35]. Both of these cases often arise in various parts of the parameter 
space of dark matter models. 
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Figure 1: Breit-Wigner enhancement [12] and the PAMELA positron excess. The analysis 
presented here is given in the two component Dirac-Majorana U{l)x and U{l)x x U{l)c 
models assuming a relic density decomposition of the Dirac and Majorana as given in the 
figure. The dominant contribution to the positron fiux comes from the annihilation of the 
Dirac particles in the galaxy. 



In the numerical analysis we fix p©^^ = 0.18 GeV/cm^ and take Pq^^ G (0.18—0.25) GeV/cm^. 
We allow the neutralino relic density to lie in the range (0.035, 0.065) and find relatively good 
fits to the WMAP and PAMELA data for both the U{l)x and the U{l)x x U{l)c model with 
total relic density in the range {Q^h'^)o + {Q^,h'^)o = (0.08, 0.12). The analysis of Fig. (II]) shows 
the PAMELA data and the positron fiux ratio in the U (l)x model and in the U {l)x x U{l)c 
model consistent with assumed densities discussed above. In this fit the dominant contribu- 
tion comes from the annihilation tptp — )■ — )■ e~^e~ and only small boost (clump) factors are 
used here, i.e., = (2 — 5). Thus, the Breit-Wigner enhancement [12] plays an important 
role in achieving a simultaneous fit to the relic density and to the positron excess. Indeed, 
the annihilation near a Breit-Wigner pole gives a significant enhancement to the annihilation 
cross section of Dirac dark matter in the galaxy obviating the necessity of using large boost 
factors. At present, we are guided by the WMAP and PAMELA data on the mass splittings 
between the Dirac component of dark matter and the vector boson mass. Generically the 
required sphtting is 2M^ — Mz' ~ (10^ — 10^) MeV depending on the leptophilic nature 
of the models and narrowness of the resonance. In Fig.([l]) the U{l)x x U{l)c has an 80 
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MeV mass splitting and the U{l)x model has 1300 MeV mass splitting. We note that since 
the U{l)x model is less leptophilic than the U{l)x x U{l)c model, we have used from the 
parameters of the model a smaller annihilation cross section for the U{l)x model than for 
the U{l)x X U{l)c model throughout the analysis as indicated by the positron ratio in the 
figure. This is motivated by the constraints from the antiproton fluxes. Thus, the p flux 
takes the form 

MT)-''i^^J.RiT)^i..),,,J^ (45) 

where T is the kinetic energy, and R{T) has been fit as in Ref. for various profile/diffusion 
models and background estimates have been obtained in [17] . The antiproton flux observed 
at the top of the atmosphere including solar modulation can be accounted for by replacing 
^p{T) — $p (T + and including a kinetic energy correction ratio. We take the Fisk 

potential ^i;' as 500 MV. In Fig. ([2]) we exhibit an analysis of the antiproton flux (signal 
plus background) in the U{l)x model. The analysis is done with the Navarro-Frenk- White 
(NFW) median model while the minimum diffusion model is unconstrained by the p data 
and is not shown. The PAMELA data exhibited in Fig. ([2]) is taken from [JB]. For the U{l)x 
model the anti-proton flux overshoots a little bit beyond E = 10 GeV but still lies within 
the limits of acceptability. For the U{l)x x f/(l)c niodel, the Z' and Z" are both leptophilic 
and there are no annihilations of ipip into qq. Because of the absence of qq final states in 
the annihilation, there is no contribution to the antiproton flux from the annihilation of the 
Dirac component of dark matter, thus the prediction of the model is not observable above the 
background. Finally, we look at the photon flux. In the angular region AQ the (differential) 
photon flux (sometime denoted d^^/dE) is given by 

For the U{l)x model there are three contributions to the photon flux. These arise from the 
qq, rf, and from bremsstrahlung (see i.e. [1511^ ). For the U{l)x x U{l)c model since the 
Z' and Z" are leptophilic with allowed final states being only in the first two generations of 
leptons, there are no final states of the type qq and rf for the Dirac component. However, 
there is an emission of continuum radiation for the Dirac component (see e.g., [50] using 
PYTHIA [Slj) which arises because ipip annihilate into ee and fifi and there is an associated 
photon continuum radiation from bremsstrahlung. 
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Figure 2: Antiproton flux in the U{l)x (signal and signal plus background). The analysis 
is done with the NFW median model while the minimum diffusion model is unconstrained 
by the p data [16] and is not shown. For the U{l)x x U{l)c, leptophilic model the p flux 
from the Dirac component of dark matter does not contribute to a signal. No boost factor 
from clumping is taken. In both models the Majorana flux is highly suppressed relative to 
the Dirac flux, and is thus not shown separately. Annihilation cross sections and local dark 
matter densities are those used in Fig.([l]). 

In Fig. ([3]) we give an analysis of the continuum photon flux in the angular region where 
the integral over the line of sight is rather insensitive to the details of the dark matter 
distribution [52] 153] (for a recent analysis with focus on the galactic center see [51]). The 
analysis is given for both U{l)x and U{l)x x U{l)c models using an isothermal proflle. The 
continuum photon flux for the U {l)x model arises mostly from qq and rf at low energies while 
the flnal state radiation, i.e. e+e~7 and /i"'"/i~7 takes over at high energies where E^/M^ — 1. 
Also shown is the EGRET [55] data and the more recent FERMI-LAT data [56] as well as 
the background flux in 10-20 region as estimated in the GALPROP analysis of [57] . For the 
U{l)x X U{l)c model the total photon flux is the suppressed contribution from the Majorana 
and the dominant Dirac source arising from the bremsstrahlung from the flnal states e"'"e~7 
and from fi~^fi~'-f. One flnds that the continuum spectrum with bremsstrahlung is in accord 
with the current experimental data. Regarding the monochromatic photon radiation from 
the annihilation of dark matter, it is suppressed by < 10""^ for the U{l)x model and the 
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Figure 3: Photon flux in the U{l)x and U{l)x x U{l)c models. The photon flux for U{l)x 
includes contributions from the quarks and taus and bremsstrahlung, while the photon flux 
for the U{l)x x U{l)c model is highly suppressed at low energies and peaks at larger energies 
from the bremsstrahlung. Annihilation cross sections and local dark matter densities are 
those as in Figs. (fT][2|) . 

prediction for this model is far below the current experimental limits. For the U{l)x x U{l)c 
model there is no coupling of the hidden sector Dirac particles to the photon if the mixing 
with the hypercharge gauge boson vanishes. Thus, at the tree level there would be no 
emission of monochromatic radiation in the annihilation of dark matter in this model. In 
the case of a small or nonvanishing mixing with the hypercharge as in the U{l)x model, this 
emission is also suppressed. 

5 CDMS-II and XENON 

The CDMS-II results mentioned in Sec.([l]) raise the possibility that 1-2 dark matter events 
may have been seen in the CDMS-II detector, and this possibility has led to a significant 
theoretical activity |58]. Many analyses within supersymmetry assume the supersymmetric 
cold dark matter is entirely composed of neutralinos. We analyze the event rates in CDMS-II 
and in XENON detectors for the case when roughly only half of the dark matter is constituted 
of neutralinos, a situation which holds for both the U{l)x and the U{l)x x U{l)c models. 
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In the analysis we use MICROMEGAS [SH] and impose the electroweak symmetry breaking 
constraints as well as all the current experimental constraints such as on g — 2, flavor changing 
neutral currents, i.e., b ^ S'-f and Bg — )■ fi'^fi" branching ratio constraints, (see i.e. [60], [61] 
and [02] for recent analyses) and require that about half the relic density as given by WMAP 
be given by neutralinos. 
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Figure 4: An analysis of the spin-independent cross section for the parameter space of 
supergravity models with nonuniversalities in the gaugino sector. The analysis given above 
is valid for both the U{l)x and U{l)x x U{l)c models with a neutralino component and 
a Dirac component of dark matter. Models are labeled by the NLSP, which under the 
constraints of radiative breaking, mass limits, and flavor changing neutral currents allow 
chargino, stau, stop, and gluino NLSPs. By far, here the chargino NLSP arises most often. 

Next we note that in the direct detection experiments, the Dirac component does not 
give an appreciable contribution and essentially the entire contribution to the event rates 
arises from the Majoranas. Speciflcally consider the U{l)x model. Here the event rates 
arising from the scattering of Dirac particles from nuclear targets are suppressed by a factor 
of relative to what one would flnd in the scattering of Majoranas. This is easily seen as 
follows: Dirac dark matter interacts with quarks in the target particles by the exchange of 
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7, Z' . The couplings of 7, Z to Dirac dark matter are suppressed by a factor of e since 
Dirac dark matter resides in the hidden sector. Thus, the cross sections arising from the 
exchange of 7, Z are suppressed by a factor of e^. Next we consider the exchange of Z' . The 
couphng of Z' [which is mostly a U{\)x gauge boson] with the Dirac component of dark 
matter is assumed to be normal size, i.e., 0{gx) ~ 0{g2)- However, its coupling with quarks 
is suppressed by a factor of e. Thus the exchange of Z' also gives a scattering cross section 
which is suppressed by O(e^). Since < 10"^ the Dirac component gives a negligible 
number of events in the direct detection experiments relative to what the neutralinos give as 
far as the Z' and Z" poles are concerned. A concrete analysis of events rates for the Dirac 
particles in these experiments has recently been given |63] and our estimate is in accord with 
this analysis. The contribution from the photon pole depends on the cutoff at small angles. 
Further, an analysis of the event rates is subject to absorption both by the atmosphere as 
well as by dirt and rock in Earth before the milli charged particle gets to the detector (see, 

e.g. m)- 

For the U{l)x x U{l)c model, the Dirac particles have no interaction with the quarks, 
so the contribution of the Dirac particles to event rates in the direct detection experiments 
is absent. Thus in either case the dominant contribution to event rates in experiments for 
the direct detection of dark matter comes from neutralinos. In Fig.dlj) we give an analysis of 
spin-independent cross section in NUSUGRA models for the parameter space of supergravity 
models with nonuniversalities in the gaugino sector so that rria = mi/2(l + ^a) with (a = 
1, 2, 3). In the analysis the NUSUGRA parameters are chosen in the following range: mo < 
3 TeV, mi/2 < 400 GeV, 6a=2,3 lie in the range (—1,1), (which statistically favors the low 
LSP mass region which is also the region of interest in this analysis) lAo/mol < 4, and 
tan /3 = (1 — 60). The current limits from CDMS, XENON, and from other experiments 
are also exhibited. We are assuming the neutralinos are contributing roughly half the relic 
abundance and roughly half the local density of dark matter. There is no rescaling by the 
dark matter density in these figures. Note the models are dominated by chargino NLSPs 
[65l [66] and the presence of a low mass chargino wall [67]. The relatively empty region in 
the range of (70-90) GeV follows from the constraint on the chargino mass being larger than 
100 GeV, and an inability for the chargino and LSP to therefore coannihilate in this region, 
along with the constraints that the stop and stau are larger in mass than 100 GeV, and the 
gluino should be larger than about 300 GeV. Further, the lightest CP even Higgs has been 
constrained to lie higher than 110 GeV. The region of low mass is mostly controlled by the 
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Table 1: Top section of the table: 


A sample set of NUSUGRA models which produce the 





Majorana component of dark matter and makes up about half the relic density of the universe 
for both the U{l)x and the U{l)x x U{l)c two component models. Middle section of the 
table: Masses for light sparticles including the neutralino Xi the light chargino x^i the gluino 
(Jj), the light stau (fi), the light stop (ti), and the CP even Higgses h,H/A (charged Higgs 
is slightly heavier) for the same set of inputs as given in the top table. Bottom section 
of the table: The Majorana relic density {Qh'^)^, {(tv)^^^^ <^ ('^'^)haio' spin-independent 
neutralino proton cross section a{SI)^p, and event rates/kg/day in germanium and xenon 
detectors corresponding to entries in the top table with an assumed 30% efficiency and with 
Pq,x ~ (l/2)p0,totai ~ 0.18 GeV/cm3. The analysis is done with a top pole mass at 171 GeV 
and the models show stability in the relic density with small changes in the pole mass. The 
(t{SI) given here are also within the range of XENON-100 and are on the edge of the limits 
recently reported in Ref. [H]. More low mass models can be seen in Fig.(jl]). 
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poles of the MSSM, while the higher mass region above 100 GeV is controlled mostly by 
coannihilations. 

6 Collider Signals 

We discuss now the collider implications of the U{l)x and the U{l)x x U{l)c models. In 
Table ([1]) we give some concrete models which generate half the relic abundance from the 
neutralino dark matter. These models produce event rates in germanium and in xenon at 
detectable levels with a relatively light spectrum. The predicted spin-independent elastic 
WIMP-nucleon cross sections are on the edge of the limits reported by XENON-100 |T3]. 
Specifically, all of the models listed in Table [1] have a light neutralino mass and several also 
have light Higgses (for recent work relating to light Higgses and the COGENT [68] and 
DAMA [69] data, see [70], [71] and [721 [73]). Further, essentially all the models in the table 
have a gluino lying in the sub-TeV range and typically all the charginos are light. However, 
some of the models have rather large scalar masses in the (1.5 - 2) TeV region indicating that 
they originate on the hyperbolic branch of radiative breaking of the electroweak symmetry 
|74j . The models listed in Table [T] share the property that the gauginos in all cases are rela- 
tively light. Thus, such models should give rise to detectable signals in the form of leptons 
and jets and missing energy at the LHC with modest luminosity (though the missing energy 
may be difficult to estimate in early runs). The models with very light gluinos could surface 
with less than 1 fb~^ at LHC center of mass energies ^/s = 7, 10 TeV, (for recent analyses 
see [75], [19] and [76]) while many of the models should be discoverable with 0(10) fb~^ 
at ^/s = 10 TeV. (In fact one can glean this from the analysis of the first listing of Ref. 
[58].) As many of the candidate models have rather light gluinos with a chargino NLSP, such 
models likely will produce missing energy which is very SM like. Large event rates can arise, 
however, from multijets, and, in particular, from b jets. One also expects a sizable amount 
of leptons in these models. For the cases where stau-coannihilation survives, the leptonic 
signals are likely to be stronger and the missing energy larger than for the chargino NLSP 
cases. However, since these models have very low SUSY scales, most of them should indeed 
be discoverable (and likely rather early) at the LHC. Additionally for the U{l)x x U{l)c 
model Z" offers the possibility of discovery at an NLC as it will have distinct signatures. Its 
decay width is significantly smaller than what a GUT-type Z' with the same mass will have. 
Additionall, it has visible decays only into ee and fifl along with radiation at the NLC. Thus, 
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a Z" of this type can be detectable at an NLC because of its distinct signatures. However, 
a full simulation of collider signals requires a separate dedicated analysis. 

7 Conclusion 

In this work we have proposed a new class of models with dark matter consisting of two, three, 
or even four components. We considered the two component model consisting of Dirac and 
Majorana particles in detail. We showed that this two component model can fit the positron 
excess seen in the PAMELA experiment as well as can produce detectable signals in the cur- 
rent direct detection experiments while satisfying WMAP relic density constraints. Thus, 
the Dirac component of the two component dark matter model allows a fit to the PAMELA 
data via annihilation of the Dirac particles close to a Breit-Wigner pole. On the other hand, 
the Majorana component of dark matter plays the dominant role in the generation of events 
in dark matter detectors. Specifically, we showed that in the two component picture it is 
possible to generate events of size 1-2 in 612 Kg-d of data in the CDMS-II detector as well 
as event rates that can be tested by the results of XENON-100, and an observable number 
of events in other ongoing direct detection experiments. Further, it was shown that models 
which lead to detectable signals in direct detection experiments are typically associated with 
a relatively light spectrum which is discoverable at the LHC with modest luminosity. Fur- 
ther, one class of models discussed in this work produces a Z" vector boson which has visible 
decays only to e+e~ and [i^ ^~ ■ The proposed models contains massive scalar fields which 
are also possible candidates for dark matter. Thus these spin zero fields in combinations 
with Dirac and Majorana particles present the possibility of a multicomponent dark matter. 
Finally, we note that it would be very interesting to investigate phenomena where both com- 
ponents play a significant role, i.e., regions of the parameter space of the model which allow 
both dark matter candidates to appear in the data analysis. Such a possibility may appear 
in certain decay fragments at the LHC where the missing energy signals from the two dark 
matter particles would be different because of their different masses and interactions. 
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